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We report the first observation of K*(892)° — > irK in relativistic heavy ion collisions. The 
transverse momentum spectrum of (K*° + K )/2 from central Au+Au collisions at y/s NN = 130 
GeV is presented. The ratios of the K*° yield derived from these data to the yields of negative 
hadrons, charged kaons, and <f> mesons have been measured in central and minimum bias collisions 
and compared with model predictions and comparable e + e~ , pp, and pp results. The data indicate 
no dramatic reduction of A'* production in relativistic heavy ion collisions despite expected losses 
due to rescattering effects. 

PACS numbers: 25.75.Dw 



Modification of meson resonance production rates and 
their in-medium properties are among the proposed sig- 
nals of a possible phase transition of nuclear matter to 
a deconfined plasma of quarks and gluons in relativistic 
heavy ion collisions Q. For resonances like K*° with 
a lifetime comparable to the time scale for evolution of 
the dense matter created in such collisions, characteris- 
tic properties such as width, branching ratio, yield, and 
transverse momentum spectra are expected to be sensi- 
tive to the dynamics and chiral properties of the high 
energy density medium which is produced ||, ^ . 

More generally, the study of short-lived hadronic res- 
onances as a means to utilize the extended spectrum 
of hadronic matter to probe hadron production under 
extreme conditions also provides important insight into 
the relative probability that a quark-antiquark pair will 
form a vector resonance meson as compared to it pseu- 
doscalar partner ||. This relates directly to the role of 
spin in hadron production in strongly interacting mat- 
ter under extreme conditions. Additionally, the study of 
higher level resonances affords a better understanding of 
feed-down to stable particles from resonance decays Q 
to further constrain thermal models of particle produc- 
tion P, pj, [t]] in nucleus-nucleus collisions. 

Resonances which decay into strongly interacting 



hadrons in the dense matter are less likely to be re- 
constructed due to rescattering of the daughter parti- 
cles. Resonances with higher px have a larger probabil- 
ity of decaying outside the system and therefore of be- 
ing detected. Alternatively, the resonance yield could be 
increased during the rescattering phase between chemi- 
cal freeze-out (vanishing inelastic collisions) and kinetic 
freeze-out (vanishing elastic collisions) [||, ^, |l0) via elas- 
tic processes like ttK — > K*° — ► irK. This regeneration 
mechanism would partially compensate for resonance de- 
cays if the expansion of the produced matter took a rela- 
tively long time (^20 fm/c), increasing the observed ratio 
of K*°/K. By systematically comparing the yields and 
transverse momentum distributions of resonances with 
other particles, it should be possible to distinguish differ- 
ent freeze-out conditions |ll , [l^] , such as sudden freeze- 
out or a slow expansion of the final state hadrons. 

Another reason the study of the K*° is interesting is its 
strange quark content. The enhancement of strangeness 
production in heavy ion collisions has long been predicted 
to be a signature of the formation of a deconfined quark- 
gluon plasma (QGP) |l3). The combined measurement of 
the K*° and (f> mesons provides an additional, unique tool 
to distinguish various hadronic expansion and freeze-out 
scenarios §, [ll], Q. 
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The if*(892)° and its antiparticle are the dominant 
resonances in the Kn system In previous rela- 

tivistic heavy ion experiments the observation of these 
resonances has been problematic due to backgrounds 
from other Kn partial waves [ 161, decays of higher mass 
resonances jlTjj elliptic flow fl8|], and detector limi- 
tations (particle misidentification, acceptance and effi- 
ciency, etc.). Due to the increased energy of the beams 
available at the Relativistic Heavy Ion Collider (RHIC), 
it was expected from simulation that the yield of these 
resonances would be sufficiently large for them to be ob- 
served using the mixed event method successfully used 
to reconstruct the 4> meson at RHIC Q. 

The detector system used for these studies was the 
Solenoidal Tracker at RHIC (STAR). The main track- 
ing device within STAR is the Time Projection Cham- 
ber(TPC) |20| which is used to provide momentum infor- 
mation and particle identification for charged particles by 
measuring their ionization energy loss (dE/dx). A mini- 
mum bias trigger was defined using coincidences between 
two Zero Degree Calorimeters (ZDC) which measured the 
spectator neutrons. A Central Trigger Barrel (CTB) con- 
structed of scintillator paddles surrounding the TPC was 
used to select small impact parameter "central" collisions 
by selecting events with high charged particle multiplcity. 

Data were taken for Au+Au collisions at y / s^" =130 
GeV. To achieve uniform acceptance in the pseudorapid- 
ity range studied plj , the collision vertex was required 
to be within ±95 cm of the midpoint of the TPC along 
the beam direction. Approximately 440,000 central and 
230,000 minimum bias events were used in this analysis. 
Particles were selected based on their momenta (p) , track 
quality, and particle identification from the TPC dE/dx. 
Since the daughters of K*° decays originate at the in- 
teraction point, tracks were selected whose distance of 
closest approach to the primary interaction vertex was 
less than 3.0 cm. Charged kaons were selected by requir- 
ing their dE/dx to be within two standard deviations 
(2er) of the expected mean. A looser dE/dx cut of 3<r 
was used for pions. Kaons and pions were required to 
have transverse momenta (pt) between 0.2 and 2 GeV/c 
to enhance track quality (21j . In addition, the daughters 
were required to have pseudorapidities |?7| < 0.8 with an 
opening angle of > 0.2 rad between them. 



The decay channels K* 



K+ and K 



K~ 



both of which have a branching ratio of 2/3, were se- 
lected for the measurements. Due to limited statistics, it 
was necessary to combine these spectra. Therefore, the 
K*° yields presented in this paper correspond to the av- 
erage value of K*° and K unless otherwise specified. 
To measure these yields, the invariant mass was calcu- 
lated for each oppositely charged Kn pair in an event. 
The invariant mass distribution derived in this manner 
was then compared to a reference distribution calculated 
using uncorrelated kaons and pions from different events. 
Typically, three or more events with similar multiplicity 
and collision vertex locations (|AZ| < 20 cm) were used 
for this "mixed-event" technique. 
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FIG. 1: a. Kn invariant mass distribution from same-event 
pairs (symbols shown every 50 MeV) and mixed-event pairs 
(histogram) from central collisions for 0.4 GeV/c < pr <2.8 
GeV/c. b. K*° invariant mass distribution after subtraction 
of the mixed-event reference distribution. A Breit-Wigner 
functional form (solid curve) was fit to the peak assuming a 
linear background residual (dashed line). The mass and width 
of the resonances used for the fit were fixed from the Particle 
Data Book j^]. The data points reflect a bin size on the 
x-axis of 10 MeV per bin. 



From the 440,000 events in the 14% most central event 
data sample, more than 1.4 x 10 10 oppositely charged 
kaon and pion pairs were analyzed. The correspond- 
ing invariant mass distribution is shown in Fig.[l].a along 
with the mixed event reference distribution. The two 
distributions were normalized to each other at Mkw — 1 
GeV/c 2 which is close to the mass region of interest for 
this measurement. The two distributions are observed to 
match well; when subtracted the resulting distribution 
exhibits a K*° signal which is approximately 15 stan- 
dard deviations above the background (FigjlJ.b). The sig- 
nal to background ratio before background subtraction 
is about 1/1000 for central events and 1/200 for mini- 
mum bias Au+Au interactions. These ratios are signifi- 
cantly smaller than the value of 1/4 observed for proton- 
proton interactions at the CERN Intersecting Storage 
Rings (ISR) p2] j, indicating the increased difficulty of 
this measurement in the high multiplicity environment 
typical of relativistic nucleus-nucleus collisions. 

As mentioned previously, higher mass resonant states 
in the Kn system as well as nonresonant Kn correlations 
also contribute to the same-event spectrum. In addition, 
particle misidentification of the decay products of the 
p, u>, r), and cause false correlations to appear in the 
same-event spectrum which are not present in the mixed- 
event spectrum used to estimate the background. Com- 
parison of the real invariant mass distribution to a ref- 
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FIG. 2: The transverse mass mr spectrum of (K*° + K* )/2 
within \y\ < 0.5 for the 14% most central Au+Au interactions 
was studied. K*° resonances having 0.4 GeV/c< pr <2.8 
GeV/c were detected. Error bars are statistical only. 



erence distribution derived using the HIJING [£3| event 
generator suggests that the residual correlation near the 
K*° mass peak may be due to the above sources. How- 
ever, accurate determination of the magnitude of this 
residual correlation requires a detailed knowledge of the 
particle production and phase space distributions for the 
above particles, including those for the p and u> which are 
presently unmeasured. Several functional forms, includ- 
ing linear and exponential, were used to fit the residual 
background in Fig.Q.b. The choice of normalization for 
the mixed-event spectrum was also varied in order to 
study the stability of the resulting K*° yield. The re- 
sulting differences in yield were within 20% in all cases, 
which was taken as a measure, in part, of the systematic 
uncertainty. 

The uncorrected number of K*° was calculated by in- 
tegrating the Breit-Wigner function fit to the data as- 
suming the linear residual background shown in Fig. |l|.b. 
In order to determine the yield, detector acceptance and 
efficiency corrections were applied as well as a correction 
for the branching ratio. This was done by embedding 

simulated kaons and pions from K*° , K decays into real 
events using GEANT, and passing them through the full 
reconstruction chain The acceptance and efficiency 
factor e depends on centrality, pr, and the rapidity of 
the parent and daughter particles. It varied from about 
10% for parent pr — GeV/c to approximately 35% for 
parent pr — 2.0 GeV/c. 

Fig. |^ shows d 2 N / (2nmTdmTdy) as function of mx — 

mo = \Jvt + m o ~ m o where m = 0.896 GeV/c 2 is 
the mass of the K*° resonance |]l7[ . An exponential 
fit was used to extract the K*° yield per unit of ra- 
pidity around mid-rapidity, as well as the inverse slope 
(T). The fit yielded dN/dy = 10.0 ± 0.9(stat) and 
T = 0.40 ± 0.02(stat) GeV for central collisions. The 
systematic uncertainty in dN/dy and T was estimated 
to be 25% and 10% respectively due to uncertainty in 
the tracking efficiency and in the determination of the 
background. Due to limited statistics, the inverse slope 



parameter derived for the central event sample was also 
used to extract dN / dy for the minimum bias sample. The 
result is 4.5 ± 0.7(stat)±1.4(sys). The additional sys- 
tematic error in this instance results from an estimate 
of the uncertainty in the inverse slope of the tot spec- 
trum. Combining all px bins, separate mass spectra of 
K*° and K were also fitted with a Breit-Wigner res- 
onant function plus a linear residual background. The 
ratio of ~K*°/K* Q = 0.92 ± 0.14(stat) was obtained for 
central events. Consequently, the average of the com- 
bined K*° and K spectra should accurately represent 
K* (892)° production within our statistics. This ratio is 
similar to K~ / K + ratio [§4). 

The K*°/h~ ratio for the top 14% most central 
collisions is 0.042 ± 0.004(stat)±0.01(sys) and 0.059 ± 
0.008(stat)±0.019(sys) for minimum bias collisions where 
h~ is the total negative hadron yield with |?7| < 0.5 
These results can be compared with K*° /h~ = 0.036 ± 
0.002 from e+e" collisions [[§ |l| || at y/s = 91 GeV and 
K*°/n- = 0.057 ± 0.009 ±0.009 from pp collisions [f§ at 
y/s = 63 GeV. The ratio K*°/h~ is observed to be ap- 
proximately constant from low to high multiplicities at 
RHIC and is compatible with that measured in elemen- 
tary particle collisions (e + e~, pp, and pp). 

With respect to studies of freeze-out conditions, the 
ratio K*°/K is more interesting and less model depen- 
dent than K*°/h~ since both particles have similar quark 
content and differ only in their spin and mass. By sim- 
ple spin statistics, the vector meson to meson (pseu- 
doscalar -(-vector) ratio would be 0.75. However, the mea- 
sured ratio is much smaller in elementary collisions [|J. 
The charged kaon results used here are an average of K + 
and K~ in the same centrality range from Ref. J24j. The 
result, K*°/K = 0.26 ± 0.03(stat)±0.07(sys) in central 
Au+Au collisions at RHIC, can be compared with the 
average value of 0.37 ± 0.02 from e + e~ & ^5), p6[ , 
pp |^7| and pp [^2[ ^8[ ^9| as shown in FigTg. 

In elementary collisions, the (p/K*® ratio measures 
the strangeness suppression to good approximation since 
there is only a small mass difference between the (j) and 
the K*°, but the strangeness quantum number for these 
particles differs by one unit, strangeness being hidden in 
the <fi meson. This ratio seems to increase in elemen- 
tary particle collisions as a function of center-of-mass 
energy(-y/s). In this study, it was found that cj)/K* a — 
0.49±0.05(stat)±0.12(sys) for the 14% most central col- 
lisions (l^] . This result is greater than that from elemen- 
tary collisions as shown in Fig. ||. The increased ratio at 
RHIC, compared to that in elementary processes, may be 
indicative of strangeness enhancement and/or additional 
effects (e.g. rescattering, coalescence |lj]) on resonances 
in heavy ion collisions. 

Since the lifetime of K*° is comparable to the time 
scale for evolution of the Au+Au collision system, the 
K*° survival probability must be accounted for when 
comparing the results from Au+Au collisions with those 
from elementary particle collisions or with thermal model 
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FIG. 3: The K M to charged kaon and <f) to K* u ratios for dif- 
ferent colliding systems as a function of i/s. Data are shown 
with quadratically combined systematical and statistical er- 
rors. The data are from collisions of e + e~ at y/s of 10.45 GeV, 
29 GeV and 91 GeV §, §, || §|, pp at 5.6 GeV @, and pp 
from the ISR §§ || at 63 and 52.5 GeV and NA27 §§ at 
28 GeV. Ratios shown are for total integrated yields except 
for the present measurements (\y\ < 0.5) and those from the 
ISR p2J at -^5=63 GeV where the ratio is for midrapidity. 



fits at chemical freeze-out. In general, the K*° survival 
probability depends on the duration of the source (At), 
the size of the source for particle emission, and the px of 
the parent K*°. If it is assumed that the difference be- 
tween the K*°/K ratio in heavy ion collisions and that 
observed in collisions of simpler systems is due to this sur- 
vival probability alone, the indication would be that At 
(<4 fm/c) is small. For large At (>20 fm/c) O without 
K*° regeneration, the K*° production should be an or- 
der of magnitude lower than the observed result jllj , and 
the low pt part of the transverse momentum distribution 
should be suppressed resulting in a larger effective inverse 
wit slope. Although the measured K*° inverse slope is 
larger than that of the charged kaons (most likely the 
result of radial flow || p4j), it is still similar to that for 
the (f> (T = 379 ± 50(stat)±45(sys) MeV) g§. This is 
consistent with the interpretation of a short time (small 
At) between chemical and kinetic freeze-out. 

Alternatively, elastic processes such as ttK — » K*° — > 
ttK are operative between chemical and kinetic freeze- 
out and partially regenerate the K*° until kinetic freeze- 
out. In a statistical model description ||, (tJ, the 
measured K*° should reflect conditions at kinetic freeze- 
out rather than at chemical freeze-out if there is a 
long lived phase in which significant rescattering takes 
place. Within the framework of this type of model, 
reasonable values || ^, |jj of the chemical and ki- 
netic freeze-out temperatures T c h and T t h, pion chem- 
ical potentials /i^ at kinetic freeze-out jl], pi, and the 
mass difference (Am = 0.4 GeV) between K*° and 
K can be obtained. These result in the ratio of 



(K*° I K) t h at kinetic freeze-out and (K*° /K) c h at chem- 
ical freeze-out being roughly (K*°/K) th /(K*°/K) ch = 
exp [(—Am + fJ-Tr)/T t h + Am/T c i] which is in the range 
of 0.3 to 1.2. The value measured in the present study 
(0.26) for kinetic freeze-out and the value of 0.37 as- 
sumed for chemical freeze-out are in the ratio of 0.7±0.2, 
which selects specific trajectories in the phase diagram 
of Hit/he and T The ratio at chemical freeze-out 
is based on a statistical model H and upon colli- 

sions of lighter systems where the chemical and kinetic 
freeze-out processes are naturally overlapped. 

Present measurements are not consistent with a long 
expansion time (Ai>20 fm/c) without significant K*° re- 
generation. They are consistent with a sudden freeze- 
out interpretation (Ai<4 fm/c). The simple estimates 
made in this paper illustrate how the study of reso- 
nance can provide important information on the dynam- 
ics and evolution of the matter produced in relativis- 
tic nucleus-nucleus collisions. More sophisticated anal- 
yses (e.g. comparison with transport models such as 
UrQMD |30| ) with improved uncertainty in the measure- 
ment of the (f> and K*° , as well as the measurement of ad- 
ditional resonances are needed to determine the evolution 
of the system resulting from central heavy ion collisions 
in detail. 

In conclusion, we have presented the first measurement 
of K* (892)° and 7?* (892)° in relativistic nucleus- nucleus 
collisions. The K*° mr spectrum from the 14% most cen- 
tral Au+Au collisions results in an inverse slope param- 
eter similar to that measured for the 4> meson in similar 
centrality. The measured yield, dN/dy = 10.0±0.9(stat) 
±2.5(sys) is relatively high compared to elementary col- 
lisions and thermal model predictions, considering the 
short K*° lifetime (cr ~ 4 fm) and expected losses due to 
rescattering of the decay daughters in the dense medium. 
The results of this study are consistent with two possi- 
ble scenarios for the dynamic evolution of the system: 
(1) a short time duration between chemical and kinetic 
freeze-out (i.e. sudden freeze-out), or (2) a long period 
of expansion characterized by high hadron density and 
significant K*° regeneration along specific trajectories in 
the phase diagram. Studies of strongly decaying resonant 
states like the K*° open a new approach to the study of 
relativistic nucleus-nucleus collisions. 
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